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Abstract 

This paper presents the Suzaku results obtained for the Sagittarius (Sgr) C region using the concept of 
X-ray reflection nebulae (XRNe) as the echo of past flares from the super massive black hole, Sgr A*. The 
Sgr C complex is composed of several molecular clouds proximately located in projected distance. The 
X-ray spectra of Sgr C were analyzed on the basis of a view that XRNe are located inside the Galactic 
center plasma X-ray emission with an oval distribution around Sgr A*. We found that the XRNe are largely 
separated in the line-of-sight position, and are associated with molecular clouds in different velocity ranges 
detected by radio observations. We also applied the same analysis to the Sgr B XRNe and completed 
a long-term light curve for Sgr A* occurring in the past. As a new finding, we determined that Sgr A* 
was experiencing periods of high luminosity already ~500 years ago, which is longer than the previously 
reported value. Our results are consistent with a scenario that Sgr A* was continuously active with sporadic 
flux variabilities of Lx = 1-3 x 10 39 erg s _1 in the past 50 to 500 years. The average past luminosity was 
approximately 4-6 orders of magnitude higher than that presently observed. In addition, two short-term 
flares of 5-10 years are found. Thus, the past X-ray flare should not be a single short-term flare, but can 
be interpreted as multiple flares superposed on a long-term high state. 
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1. Introduction 

In the Galactic center (GC), prominent 6.4 keV line of 
neutral iron (Fe I Ka) has been detected from giant molec- 
ular clouds (MCs) such as Sagittarius (Sgr) B (I ~ 0.°7; 
Koyama et al. 1996), Sgr C (I ~ 359.°5; Murakami et al. 
2001; Nakajima et al. 2009), and Sgr A (I ~ 0.°1; Park et 
al. 2004; Muno et al. 2007). The spectra show large equiv- 
alent width (EWqa keV > 1 keV) and strong absorption 
(N-a > 10 23 H cm~ 2 ) to the continuum, which suggests 
that the 6.4 keV bright MCs are due to the irradiation 
and fluorescence by possible external X-rays (X-ray re- 
flection nebulae: XRNe) rather than particle irradiation 
(e.g., Yusef-Zadeh et al. 2007). Subsequently, the XRN 
scenario has become more conclusive owing to the discov- 
ery of the short-term (a few years) 6.4 keV variability in 
the small (a few light-years) regions in Sgr B (Koyama 
et al. 2008) and Sgr A (Ponti et al. 2010). Moreover, in 
Sgr B, the same time variability of hard continuum X-rays 
(E > 8 keV) has been found to be in correlation with the 
6.4 keV line by Terrier et al. (2010) and Nobukawa et al. 
(2011). These results promoted the idea of the external 
irradiation source being a past flare of the super massive 
black hole Sgr A* (e.g., Koyama et al. 1996; Sunyaev et al. 
1993; Murakami et al. 2001; Ponti et al. 2010; Nobukawa 



et al. 2011). 

In the XRN context, one can derive the X-ray light 
curve of Sgr A* for the past several hundred years. The 
flare luminosity and look-back (delay) time depend on the 
position of the XRN relative to Sgr A* and the sun (ob- 
server). However, most of the previous estimates used 
solely the projected distance or indirect information of 
line-of-sight positions. 

Ryu et al. (2009) successfully developed an original 
method of determining the line-of-sight position for the 
Sgr B XRNe. Their method involves careful spectral anal- 
ysis of XRNe in combination with the spatial and spectral 
information of the GC plasma X-ray emission (GCPE). 
Uchiyama et al. (2011) recently performed an intensive in- 
vestigation on the 6.7 keV (Fe XXV Ka) line profile for the 
entire GC region (~ 5° x 2°); thereafter, Uchiyama et al. 
(2012) extended the work to other emission lines of highly 
ionized sulfur (S), argon (Ar), and calcium (Ca). These 
studies quantitatively constructed the GCPE as a two- 
temperature plasma and determined its three-dimensional 
spatial distribution. 

On the basis of these new pictures of the GCPE, we es- 
tablish the method to measure the three-dimensional po- 
sitions of the Sgr C XRNe. The observation and analyses 
of the Sgr C data are given in section 2 and 3. We derive 
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a face-on view of XRNe in section 4.1 and compare its re- 
liability with that of radio observations in section 4.2. We 
then combine these results with the Sgr B data and re- 
sults (Ryu et al. 2009) to calculate the time delay and the 
corresponding X-ray flare luminosity to reconstruct the 
long-term (^500 years) light curve for Sgr A*, which oc- 
curred in the past. The derivation method and discussion 
for the activity history of Sgr A* are given in section 4.4. 

Throughout this paper, we adopt 8.0 kpc (Ghez et al. 
2008) as the distance between the sun and Sgr A*; thus, 

1. °0 corresponds to 140 pc at the GC. The parameter un- 
certainties are quoted at the 90% statistical confidence- 
level (1.64 a) range unless noted otherwise. 

2. Observation 

We have made two deep pointing observations on the 
Sgr C region separated by ~5 years between February 
2006 and September 2010. These observation are con- 
ducted using the X-ray Imaging Spectrometer (XIS) at the 
focal planes of the X-ray Telescope (XRT) on board the 
Suzaku satellite. The XIS system consists of three sets of 
front-illuminated (FI) charge coupled device (CCD) cam- 
eras (XIS 0, 2, and 3) and one set of back-illuminated (BI) 
CCD camera (XIS 1); each CCD chip contains 1024x1024 
pixels (1 pixel = 24 /im x 24 fj,m) for a 17.'8xl7.'8 field 
of view. Two calibration sources of 55 Fe are installed to 
illuminate two corners of each CCD for absolute gain tun- 
ing. We operated the XIS in the normal clocking mode 
with a read-out cycle of 8 s. The details of Suzaku, XIS, 
and XRT are given in Mitsuda et al. (2007), Koyama et 
al. (2007a), and Serlemitsos et al. (2007), respectively.The 
observation log of Sgr C is shown in table 1. Archive data 
information on the Sgr B and other relevant GC regions 
can be found in Ryu et al. (2009) and Uchiyama et al. 
(2011). 

3. Data Reduction and Analyses 

We perform the data reduction and analysis using 
HE ADAS software version 6.12. The calibration database 
used to calculate the response of XIS and the effective 
area of XRT is the version released 1 on 2011-02-10. Since 
the relative gains and response functions of the FI CCDs 
are essentially the same, the FI spectra are merged 2 in the 
spectral fitting. For the images described in the following 
sections and spectral analyses, the non- X-ray background 
(NXB; Tawa et al. 2008) is generated using the night- 
Earth observation data and is subtracted from the raw 
data. The XRT vignetting effects and the exposure-time 
differences have been corrected. 

3.1. Images 

We construct the XIS images of the Sgr C region with 
the 6.4 keV-line (Fel), the 6.7 keV-line (FeXXV), and the 

1 http:/ /www. astro, isas.ac.jp/suzaku/caldb/ 

2 The XIS 2 suddenly became unusable on November 9, 2006, pos- 
sibly due to a micro-meteoroid impact on the CCD. Therefore, 
after this epoch, the merged FI CCD were XIS and 3. 



2-3 keV bands (see figure 1). The images are generated 
by adding data from multiple observations together (c.f., 
table 1 and Uchiyama et al. 2011) to cover a wide region 
near Sgr C. In addition, the FI and BI data are combined 
to increase the statistics. For the 6.4 keV and 6.7 keV 
line images, the fluxes underlying the continuum are es- 
timated from the 5-6 keV band and subtracted from the 
6.3-6.5 and 6.6-6.8 keV bands, respectively. Details of 
the subtraction procedures are explained in Nobukawa et 
al. (2010). 

Figure la displays many bright clumps in the 6.4 keV 
band near Sgr C. These are spatially associated with giant 
MCs (see figure lc) as observed in the radio CS (J=l— >0) 
line by Tsuboi et al. (1999). We designate the brightest 
three 6.4 keV clumps: CI, C2, and C3 (solid circles; fig- 
ure la). The clumps CI and C2 were previously identified 
as XRNe by Nakajima et al. (2009), while C3 is a new 
XRN candidate revealed in the 2010 observation. 

Figure lb (the 6.7 keV line image) traces the high- 
temperature plasma (HP) of the GCPE (section 3.2), and 
shows a more uniform distribution than the 6.4 keV emis- 
sion shown in figure la. Figure Id shows the 2-3 keV band 
image including the 2.45 keV line (Sxv) and displays the 
distribution of the low-temperature plasma (LP) and in- 
dividual thermal diffuse sources, mainly supernova rem- 
nants candidates, with temperatures of kT ~d keV. The 
faint diffuse sources near the Sgr C clumps are the super- 
nova remnants G359.41— 0.12 and its outflow "Chimney" 
(Tsuru et al. 2009). To investigate the nature of the 
GCPE (HP and LP), we selected two regions with rel- 
atively weak 6.4 keV emission for spectral analyses, and 
designate them as refl and ref2 (solid circles; figure la, b, 
and d). 

3.2. Spectral Model 

The X-ray spectra of the GC region have been exten- 
sively studied by Ryu et al. (2009) and Uchiyama et al. 
(2012). The spectra can be reproduced by the superposi- 
tion of four components: the GCPE, the XRN emissions 
(XRNE), the foreground emission (FE), and the cosmic 
X-ray background (CXB). These components have the 
following properties. The GCPE has two-temperature 
plasmas with kT— 1 keV and kT— 7 keV for LP and 
HP, respectively The XRNE is expressed as neutral ion 
lines (Ka and K/3) associated with a power-law compo- 
nent. Since the XRNE is formed by Thomson scatter- 
ing of a continuum component and fluorescence of iron 
atoms in the MCs, it shows a large absorption (Absl) of 
Ah ~ 10 23 H cm -2 (e.g., Murakami et al. 2001), which 
is approximately equal to the Ah through the MC. The 
GCPE including the HP and LP is extended in the GC 
region, and the MC is inside this region. Thus, a fraction 
(R) of the X-ray from the GCPE is not absorbed by the 
MC while the other fraction (1— R) is absorbed by the 
MC (Absl). Both GCPE and XRNE exhibit an common 
interstellar absorption (Abs2) of Ah ~ 6 xlO 22 H cm -2 
between the GC and the sun (Sakano et al. 2002; Ryu et 
al. 2009). The FE is approximated to a thermal plasma 
of kT-1 keV with small absorption. The CXB (e.g., 
Kushino et al. 2002) has the absorption of Abs2 by twice 
(interstellar absorption in front of and behind the GC), 
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Table 1. Suzaku Observations at the Sgr C Region. 



Target 


Obs. ID 


Obs. point (FOV center) 


Obs. date 


Effective 


name 




a (J2000) 5 (J2000) 


(yy-mm-dd) 


exposure 


SgrC 


500018010 


17 h 44 ra 37?30 -29°28'10."2 


2006-02-20 


106.9 ks 


SgrC 


505031010 


17 h 44 m 58f01 -29°22'51"2 


2010-09-25 


82.6 ks 




Galactic; longitude I [dcgreel 



Fig. 1. Added XIS (FI+BI) images near the Sgr C region. The images are binned with 12 X 12 pixels and smoothened with a Gaussian 
kernel of 5". (a) 6.4 keV line (Fel) image obtained by subtracting the continuum component (see text) from the 6.3—6.5 keV band 
flux. The XIS field of views (FOV; 17.'8xl7.'8) of the 2009 and 2006 observations are represented by the large boxes. The regions for 
spectral analyses (CI— 3, ref 1—2) are indicated by solid circles, (b) Same as panel (a), but for the 6.7 keV line (Fexxv) image, (c) 
6.4 keV line image (gray) with the superposed radio CS (J=l— >0) contours (Tsuboi et al. 1999) in the velocity range of —100 ±20 
km s" 1 (red), — 60±20 km s" 1 (green), and +50±30 km s _1 (blue). The contour level begins at the intensity of 5 K km s _1 and 
increases every 5 K km s — 1 . (d) Same as panel (a), but for the 2-3 keV band image (color). The region of G359.41— 0.12 in Tsuru 
ct al. (2009) is represented by the dashed circle. 



in addition to Absl (MC). Thus, the spectral model can 
be expressed by equation (1). A schematic picture of this 
model is shown in figure 2. 

f(E) = [Absl x Abs2 x (1 - R) + Abs2 x R] x (HP + LP) 
+[Absl x Abs2] x XRNE 

+[Absl x Abs2 x Abs2] x CXB + FE (1) 

Ryu et al. (2009) successfully applied this model to the 
X-ray spectra of the 6.4 keV clumps near Sgr B. In ad- 
dition, Muno et al. (2004) and Nobukawa et al. (2010) 
reproduced the X-ray spectra of Sgr A regions with sim- 
ilar models. Therefore, we apply the same model to the 
Sgr C region. In figure 3, we show a simulated spectrum 
(XIS/FI; 0.5-10.0 keV) with the typical GC parameters 
(see caption). The plasma emission code used in this pa- 



per is APEC 3 and the absorption is evaluated in solar 
abundances 4 . 

3.3. Spectral Fitting 

As shown figure 4, we extract seven spectra from five 
regions (solid circles in figure 1). The spectra of CI and 
C2 are obtained twice in the observations of 2006 and 
2010, whereas C3 is only observed in 2010. All spec- 
tra are simultaneously fitted in the band of 0.5-10.0 keV 
with the model given in equation (1). Since the spectral 
model is highly complicated with many physical parame- 
ters, it is impractical to fit the data with all parameters 

3 Astrophysical Plasma Emission Code: A model of an emission 
spectrum from diffuse gas in collisionally-ionizcd equilibrium 
(Smith ct al. 2001). 

4 The solar abundance in this paper is referred to the values from 
Anders & Grevesse (1989). 
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Fig. 2. A schematic face-on view for linc-of-sight alignment. 
Coordinates X, Z, and Y correspond to directions of I, ft, and 
line of sight, with respect to Sgr A*. 
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Fig. 3. Simulated XIS/FI spectrum with typical GC param- 
eters reported by Ryu ct al. (2009). Model components 
of GCPE (kT LP =1 keV, kT HP =6.5 keV, a=0.27, R=0.5), 
XRNE (EW 6A kcV =1.6 keV, V=1.7), FE, and CXB are 
shown in red, blue, green, and cyan lines, respectively. The 
column densities of Absl and Abs2 are set to lOx 10 22 H cm -2 
and 6xl0 22 H cm -2 , respectively. 

free. Therefore, we adopt several reasonable constraints 
in the fitting. The parameter settings are the same as 
those of Ryu et al. (2009). The FE is fixed to an APEC 
model with kT=0.85 keV and abundance of 2=0. 011, ap- 
plied by small absorption of Ah=0.17x 10 22 H cm~ 2 (Ryu 
et al. 2009). The center energy of the Fe Ka and K(3 lines 
are respectively fixed to 6.4 keV and 7.05 keV with the flux 
ratio (Ka/K/3) of 0.125 (Kaastra & Mewe 1993). The tem- 



perature of the HP (kTnp) is fixed to 6.5 keV (Koyama et 
al. 2007b), and the abundances of HP and LP are fixed to 
one solar (Z = 1). The parameter a, normalization ratio 
of HP/LP, can be approximated to a constant near Sgr C. 
This estimation is inferred from results of Uchiyama et 
al. (2012) (c.f., figure 3 and table 2 therein), which indi- 
cate that the two plasmas exhibit nearly the same spatial 
distribution with an e-folding length of ^0.°6. 

After all, the temperature of LP (kTLp), the normal- 
ization ratio of HP/LP (a= normup/normLp), the pho- 
ton index (r) of the power law, and the equivalent width 
(EWqa kov) of the 6.4 keV line are free parameters, but 
are linked in all regions including Cl-3 and refl-2. These 
settings are designed to verify a common and consistent 
XRN model. Other parameters, such as line intensity 
(^6.4 kev) and absorption (Ah), are free and independent 
for each region. To investigate the time variability for CI 
and C2, the 6.4 keV intensities of the 2006 and 2010 spec- 
tra are set as free parameters, while other parameters are 
common between the two spectra. 

As is seen in figure lb, CI is entirely overlapped with the 
" Chimney" , which is a thermal plasma outflowing from 
G359.41— 0.12. In order to include the Chimney contam- 
ination, we add the spectral model with fixed parameters 
(kT, A H , Z) reported in Tsuru et al. (2009). The flux 
(normalization of Chimney) is calculated from the average 
surface brightness and the region size of CI (radius=2f0), 
which is 35% of the GCPE flux. 

All seven spectra are generally reproduced by the 
model of equation (1) with a reduced x 2 /d.o.f of 
3773/3031=1.26. The fitting results are shown in figure 4 
and the best- fit parameters are listed in table 2. The 
large EWqa keV of ~ 1.2 keV and the strong absorption 
(A H (Absl)^ 10 23 H cm -2 ) agree with the XRN scenario 
for all the regions. The Ah of Abs2 are in the range of 
5.0-6.5 xlO 22 H cm~ 2 . These values are consistent with 
the interstellar medium (ISM) absorption toward the GC 
region (Ryu et al. 2009; Nobukawa et al. 2010; Sakano 
et al. 2002). As for the time variability of the 6.4 keV 
line between the 2006 and 2010 observations, CI exhibits 
a small increase of 8% at 2.9-cr level, while C2 shows no 
significant change. 

For confirmation, we also apply independent fitting with 
free EWqa keV parameters for each region. The spectra 
of CI, C2, and C3 show EW 6A kcV of 1.1-1.5 keV, 1.1- 
1.6 keV, and 0.7-1.3 keV, respectively. These values are 
consistent with the fitting result when EWqa kcV is a com- 
mon free parameter {EWqa k c v=l-15-1.27 keV; table 2), 
and agree with the XRN scenario of large EWqa kcV- 
However, the independent EWqa kcV fittings result in 
large statistical errors in the R parameters (about 2-4 
times larger than those in table 2), which is not practical 
for determining cloud positions (section 4.1) and follow- 
ing discussions. Thus, we adopt the fitting results with 
EWqa kov as a common free parameter between the re- 
gions (table 2). 
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Fig. 4. Fitting results for all selected spectra (regions represented by black solid circles in figure fa) near Sgr C. The Ft (XfS 0+2+3 
for 2006; X1S 0+3 for 20f0) spectra with 1-a error bars and the best-fit models (c.f. , figure 3) are represented with solid crosses and 
dashed lines, respectively. The best-fit parameters are summarized in table 2. 



4. Results and Discussions 

4-.1. Line- of -sight Position of the XRNe 

As illustrated in figure 2, the parameter R repre- 
sents the line-of-sight position of the MC in the GCPE. 



Depending on the distance from Sgr A*, the GCPE is 
more or less contaminated by the Galactic ridge plasma 
emission (GRPE), which has a similar two-temperature 
spectral structure but a more extended spatial distribu- 
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Table 2. Best-fit spectral parameters of the Sgr C regions* 









Local parameters 






Region 


Obs.Yr 


Nl (Absl) 

rl \ / 


Nl (Abs2) 


^6.4 kcV"^ 


FE norm§ 


GCPE noririLpH 


Fraction R# 


CI 


2006 
2010 


6-5±8;| 


c 1+0.3 
°- x -0.3 


-, yo + 0.05 

L - 1 °-0.05 


060+ 004 
u.uou_ 004 


O Q+0.1 


n Q q +0.03 

u - z0 -0.04 


C2 


2006 
2010 


ii-4li: 6 4 


6 5+ 03 

D.O_ 2 


i n +0.03 
1 - 01 -0.03 

1 r 1+0.03 

i - oi -0.03 


o-ioo±8:8gg 


2 0+ 01 
z - u -o.i 


72+ - 12 

u - ' z -0.03 


C3 


2010 


8 7 +0A 


c (1+0.3 
O.U_ 2 


1 1fi+ 005 


0.082i° : §°| 


O Q + 0.1 

°- u -0.1 


2fi +0 - 03 

u - zo -0.03 


rcfl 


2006 




c Q+0.2 

°-°-o.i 


u - 3 °-0.04 


0.059i°;fe 


9 q+0.1 

z -°-o.i 


40+ 004 

u -^ u -0.03 


ref2 


2010 


7.6i°; 5 7 


6 4 +0A 


n 70+0. 05 


0.08218:88" 


q q+0.1 
°-°-0.1 


31 + 05 
u - oi -0.05 



Global parameters 

kT LP a ratio (^HP) SW 6 .4keV Photon index T 

0.9218:81 keV 0.28 ±g;g| 1.19lg;gS keV 1.63l8;gg 

x7d.o.f ** = 3773/3031 = 1.26 



* The spectral regions are shown in figure 1. The uncertainties arc at 90% confidence level (statistical), 
t The value of column density in the unit of 10 22 Hem" 2 , 
t Intensity in the unit of 10 -6 photons era" 2 s _1 arcmin~ 2 . 

§ The normalization factor of the APEC model for the foreground emission. It is normalized with the region size A [arcmin 
and expressed as 7.07 X 10 -13 /(47rD 2 A) EM [cm - 5 arcmin -2 ], where D and EM are the distance to the source [cm], and the 
emission measure [cm -3 ], respectively. 

I In the same expression as §, but the normalization factor of LP in the GCPE. 

* R is the fraction of the GCPE not absorbed by MC (see equation (1) and figure 2), which is the indicator of the line-of-sight 
position for XRNe (see equation (2)). 

** The results are obtained by the simultaneous fitting of the 14 spectra from the FI and BI data. 



tion (Uchiyama et al. 2012). To estimate the MC posi- 
tions numerically, we use the spatial distributions of these 
plasmas. 

We introduce new coordinates (A, Y) with origins at 
Sgr A* (I = -0.°056,6= -0.°046; Yusef-Zadeh et al. 1999), 
in which X and Y correspond to the directions of / and 
line of sight, respectively (see figure 2). The offset of 
b is ignored because the relevant XRNe are all on or 
near the Galactic plane (i.e., \b\ <0.°1). Uchiyama et 
al. (2011) and Uchiyama et al. (2012) derived the X- 
axis flux distribution for many emission lines of ionized 
atoms (e.g., Fe and S) in addition to the soft and hard 
X-ray continuum bands. The data were fitted with a two- 
exponential function along the Galactic plane (X-axis): 
Jcexp(— X/Iq) + iR.exp(— X/I-r). The values of e-folding 
length for both HP and LP are similar in the GCPE and 
GRPE. In this paper, we adopted ?c=0.°6, Zr = 50°, and 
^c/^R=10 as the mean values (see table 2 in Uchiyama et 
al. 2012). 

We assume that the X-Y plane distribution of the plas- 
mas is an oblate spheroid expanded with two-exponential 
function as I(X,Y) = A x exp(— r/rcc) + exp(— r/rcR,), 
where r = ^X 2 + Y 2 . A is the flux ratio (GCPE/GRPE), 
and rcc and tgr are the e-folding scale lengths of GCPE 
and GRPE, respectively. Integrating I(X,Y) along the Y- 
axis as I(X) = Jl° I(X,Y) dY, we find that the resultant 
projection of A-distribution is nearly the same exponen- 
tial shape as observed by Uchiyama et al. (2012). We then 
compare I(X) with the results of Uchiyama et al. (2012) 
to match the shapes. The parameters are nicely deter- 
mined to be rcc=55 pc (or 0.°4), tgr=5000 pc (or 36°), 
and A=900. 



Considering that a MC/XRN is located at (Xq,Yq) and 
Ao as its I offset from Sgr A*, the line-of-sight position 
(Yq\ also c.f., figure 2) is determined from the following 
relation, using the obtained R value (table 2). 



R(y) f\ pc I(Xo,Y)dY 



(2) 



We convert the unit of (A, Y) into the actual length of 
light year (ly). Figure 5 illustrates the resultant face-on 
view of XRNe around Sgr A* on the X-Y plane. Although 
CI and C2 are proximately (~40 ly) located in the pro- 
jection (figure la), CI is in front of the A-axis, while C2 
is largely behind CI by ~400 ly beyond the X-axis. 

4-2. Association with the Radio Molecular Clouds 

We have shown that CI, C2, and C3 are largely sepa- 
rated in the line-of-sight positions, although they are prox- 
imately located at the projected distance. We search for 
the radio counterparts of CI, C2 and C3 from the CS 
(J=l— >0) maps (Tsuboi et al. 1999) in different veloc- 
ity ranges from —250 km s _1 to 250 km s _1 . As shown 
in figure lc, C2 and C3 are respectively in good coinci- 
dence with MCs in the velocity ranges of — 60 ± 20 km s _1 
and +50±30 km s^ 1 . On the other hand, CI is observed 
in the two separated velocity ranges of — 60 ± 20 km s _1 
and —100 ±20 km s _1 . Using the data corresponding 
C2 and C3, we derived a conversion factor between the 
X-ray column density (Absl; table 2) and the integrated 
CS intensity (contours; figure lc), which is [CS/Ah] — 
2.5 [K km s-y(10 22 H cm" 2 )]. The X-ray A H of MC 
(Absl) used in this study is determined directly from the 
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Fig. 5. Face-on view of Sgr C XRNe (filled circles) around 
Sgr A*. The parabolas (dashed lines) represent the equal-time 
delay (T,) contours for the X-ray echoes from Sgr A*. Sgr B 
results from Ryu et al. (2009) (open circles) are also added 
(see section 4.3). Parameter details are given in table 3. 

absorption of GCPE continuum behind the cloud (c.f., 
figure 2). This method is available even if the MC is not 
an XRN, and thus does not depend on the properties of 
the illuminating source. Assuming the same [CS/Ah] fac- 
tor near Sgr C, the observed X-ray column density of CI 
(6.0-6.8xl0 22 H cm -2 ; table 2) indicates a CS intensity of 
15-17 K km s _1 . This generally excludes the —60 km s _1 
MC as the counterpart of CI, because the average CS in- 
tensity is only 5-10 K km s _1 . The -100 km s _1 MC 
has a CS intensity of 10-25 K km s _1 , which is in agree- 
ment with the iVu estimated by the X-ray data. Thus, 
CI is most likely to be situated in the velocity range of 
-100 ±20 km s" 1 . 

Sofue (1995) identified two large l-V structures known 
as "Arm I" stretching from (359.°3, -150 km s" 1 ) to 
(0.°9, +80 km s" 1 ) and "Arm II" stretching from (359.°4, 
-80 km s" 1 ) to (0.°1, +60 km s" 1 ). According to the ve- 
locities, the Sgr CI XRN may belong to Arm I, while C2 is 
in Arm II. Assuming the uniform circular rotation around 
Sgr A*, Sofue (1995) predicted that Arm I and Arm II 
are respectively located in the foreground and background 
with respect to Sgr A* (c.f., figure 10 of Sofue 1995). The 
X-ray face-on view (figure 5) is basically consistent with 
the radio result. 

4-3. XRN Parameters and Two-temperature Structure of 
the Galactic Center Plasma 

The photon index and EWqa koV for the 6.4 keV clumps 
in Sgr C are 1.54-1.71 and 1.15-1.27 keV, respectively (ta- 
ble 2) . These values are approximately equal to the Sgr B 
results (Ryu et al. 2009). The large EW is in favor of the 
X-ray reflection and fluorescent origin (e.g., Koyama et 



al. 1996; Murakami et al. 2001). The photon index is con- 
sistent with the canonical active galactic nucleus (AGN; 
e.g., Ishisaki et al. 1996), which favors the XRN scenario 
owing to the past flare of Sgr A*. 

The best-fit temperature of LP (kTj_,p) and the mixing 
ratio (a = normup/norniLp) in the GCPE at Sgr C are 
0.87-0.95 keV, and 0.23-0.29, respectively. These values 
are nearly equal to those in Sgr B region: kXLp =0.81- 
0.91 keV and a =0.26-0.28 (Ryu et al. 2009). Sgr C and 
Sgr B are located at nearly the symmetrical positions with 
respect to Sgr A* on the Galactic plane, and hence these 
results indicate that the GCPE distribution is symmet- 
rical not only in the flux profile (Uchiyama et al. 2012) 
but also in the spectral shape. These results lead us to 
re-estimate the line-of-sight positions of Sgr B MCs using 
the same methods for Sgr C (section 4.1). The line-of- 
sight positions derived from parameter R of the Sgr B 
MCs (Ryu et al. 2009) are also plotted in figure 5. 

4-4- Activity History of Sgr A* 

On the basis of the XRN scenario of the past flares of 
Sgr A*, we construct the parabolas, with their common 
foci at Sgr A*, to illustrate the equi-delay time (T*) con- 
tours of the X-ray echoes (see figure 5). The time delay 
of the echo is given as below. 



r, = (v/A 2 +r 2 +y)/c [yr] (3) 

Using the fluorescent Fel intensity (Iqa kev), the dis- 
tance (D*=v / A 2 + Y 2 ), and the column density of XRN 
(Absl; table 2), the required X-ray (2-10 keV) luminosity 
of Sgr A* (L*) can be expressed as below (e.g., Nobukawa 
et al. 2008; Sunyaev & Churazov 1998). 



-v 1 

-2\-l 



= 1.0 x (Tg.4 k c v /10 photons cm 
x(L>* /300 ly) 2 x (N B /TO 22 H cm" 2 ) 

[10 40 



'cmin 2 ) 



ergs (4) 

Here, the isotropic radiation of Sgr A* with a photon in- 
dex of T = 1.6 is assumed (table 2). We then derive T* 
and for three Sgr C XRNe (Cl-3) and the Sgr B XRNe 
(Ryu et al. 2009). The results are summarized in table 3 
and the long-term light curve of Sgr A* is plotted in fig- 
ure 6. Sgr A* was active with short flares and possibly in 
a continuous high-luminosity state with = 1-3 x 10 39 
erg s -1 from the past ~50 to ^500 years ago. Thanks to 
the measurement of line-of-sight position of C2 (see fig- 
ure 5), we find that Sgr A* was experiencing periods of 
high luminosity already 500 years ago, which extends the 
light curve of Inui et al. (2009) and Ponti et al. (2010) 
back by 200-400 years. 

In the past high state, the luminosity shows sporadic 
variabilities (short flares) at ^200 and ~100 years ago as 
suggested by Sgr B (M 0.74) and Sgr CI, respectively. The 
short flare at 200-300 years ago would be responsible for 
the short-term time variability of Sgr B2 (Nobukawa et al. 
2011; Koyama et al. 2008), and hence the variability time 
scale would be 5-10 years. The other short flare at ~100 
years ago has similar time-scale as inferred by the 6.4 keV 
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Table 3. Three-dimensional parameters of XRNe and past luminosity of Sgr A*. 



XRN 


Obs. 


xt 


yt 


dI 


tJ 


Lit 


ID 


Year 


[iy] 


[iy] 


[iy] 


[year] 


[10 39 erg s" 1 ] 



Sgr CI 2010 -230 -241^0 333+^ 92+^ 3.1+Q5 

SgrC2 2010 -216 189±^° 287±^° 4761^° 1.2±^ 

Sgr C3 2010 -162 -173±fg 237±2° 64 ±io - 8 -o'i 

Sgr B2 2005 327 -55i^ 331±f 28ll^ 1.7±g;i 

Sgr Bl 2006 258 -75±|| 269±?S 198± 4 £ 1.0±g;f 

MO.74-0.09 2005 363 -202±|g 416^ 219±i$ 2.7i£ :4 

M0.74-sub 2007 363 -ggOjJ) 459 -4Q 182 -2Q 2 ;%-ol 

* The uncertainties of Y are at 90% confidence level (statistical). Errors of D t , T», and 
L„ are estimated using the Y ranges. 

t Coordinates and Y correspond to the positions in the Galactic longitude (£) and the 

line of sight from Sgr A*, respectively. 

§ Distance from Sgr A* calculated from Dt,=\/X 2 + Y' 2 . 

II Time delay of echo from Sgr A* calculated from equation (3) and counted from 2010. 

* Required past X-ray luminosity of Sgr A* calculated from equation (4), also sec text. 
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Fig. 6. X-ray light curve of Sgr A* in the past 500 years. 
Data points of Sgr C and Sgr B are shown by filled and open 
circles, respectively. Error bars contain distance uncertainties 
(see table 3). The inset panel shows the short-term variability 
(l-cr error) indicated by the two observations of CI on the 
6.4 keV line (see table 2), in which the distance uncertainties 
are not included. The present luminosity of Sgr A* is quoted 
from Baganoff et al. (2001) and Porquet et al. (2003). 

time variability (increase) of CI in 4 years (2006-2010) at 
3-er level (table 2). This flare may correspond to the time 
variability of the Sgr A XRNe. Muno et al. (2007) and 
Koyama et al. (2009) reported variability of the 6.4 keV 
line from the Sgr A cloud in a short time-scale of 3-5 years 
at a ~ 5cr sigma level. Then Ponti et al. (2010) and Capelli 
et al. (2012) predicted the flare-like events occurred at 
~ 100-400 years ago, which is basically consistent with our 
results. Ponti et al. (2010) proposed that the flare should 
have started ~400 years ago, while Capelli et al. (2012) 
suggested that a high period ended ^150 years ago from 
the illumination pattern in the Sgr A complex. However, 



the estimated flare luminosities are different: ~ 10 39 
erg s _1 was suggested by Ponti et al. (2010), while L* ~ 
10 38 erg s _1 was suggested by Capelli et al. (2012). These 
inconsistencies are mainly attributable to unclear line-of- 
sight positions and partly due to the difficult estimation 
of ./Vh for the Sgr A XRNe. The luminosity determined 
by Capelli et al. (2012) is - 10 37 -10 38 erg s" 1 around 70 
and 130 years ago, which is about 10 times lower than 
that obtained here: L* ~ 10 39 erg s _1 at ~90 years ago. 
The high luminosity indicated by CI may be attributed 
to a short-term flare that occurred in the fading period 
reported by Capelli et al. (2012). A unified analysis for 
the Sgr A XRNe as performed in Sgr C and Sgr B in this 
paper would slove these problems. 

The average X-ray luminosity in the high state was 
Lf c ~ 2xl0 39 erg s" 1 , which is only - 4 x 10~ 6 of the 
Sgr A* Eddington limit (L E ^ 5xl0 44 erg s _1 ) with a 
mass of 4xl0 6 M (Ghez et al. 2008). Still, from this 
luminosity, we may state that Sgr A* has been a low- 
luminosity AGN (e.g., M 81; Ishisaki et al. 1996) for the 
past 500 years. 

The present X-ray luminosity of Sgr A* is ~ 



10 



33-35 



erg s (Baganoff et al. 2001; Porquet et al. 



2003), and hence a sudden luminosity drop by 4-6 orders 
of magnitude should have occurred within ~100 years. 
Schawinski et al. (2010) reported that the quasar IC 2497 
with mass of ~ 10 9 Mq experienced an dramatic luminos- 
ity drop by more than 4 orders of magnitude in 45000- 
70000 years. On the other hand, for the Galactic X- 
ray binaries with a black hole mass of ~40 M Q (e.g., 
GRS 1915+10; Greiner et al. 1996), a state change from 
high/soft to low/hard was detected on a time scale of 
1 hour. Assuming that the state-change time scale is pro- 
portional to black-hole mass, at a mass of 4xl0 6 M Q , 
70000 years (for 10 9 Mq) and 1 hour (for 10 Mq) corre- 
spond to 280 years and 46 years, respectively Then, the 
time scale of the luminosity drop for Sgr A* is estimated 
to be as fast as 50-300 years, which is consistent with the 
prediction in our light curve (figure 6). 
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